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ABSTRACT 

We study the nonlinear growth of kinetic gyroresonance instability of cosmic rays (CRs) induced by large 
scale compressible turbulence. This feedback of cosmic rays on turbulence was shown to induce an important 
scattering mechanism in addition to direct interaction with the compressible turbulence. The linear growth is 
bound to saturate due to the wave-particle interactions. By balancing increase of CR anisotropy via the large 
scale compression and its decrease via the wave-particle scattering, we find the steady state solutions. The 
nonlinear suppression due to the wave-particle scattering limit the energy range of CRs that can excite the 
instabilities and be scattered by the induced slab waves. The direct interaction with large scale compressible 
modes still appears to be the dominant mechanism for isotropization of high energy cosmic rays (> 100 GeV). 

Subject headings: instabilities-MHD-scattering-turbulence-cosmic rays. 
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1. INTRODUCTION 

Cosmic rays and turbulence are essential components in as- 
trophysical systems. Their interactions are crucial for high en- 
ergy phenomena and dynamics of various systems. The reso- 
nant and non-resonant (e.g., transient time damping, or TTD) 
interactions of cosmic rays with MHD turbulence are the ac- 
cepted p rincipal mechanis m to scatter and isotropize cosmic 
rays (see ISchlick eiser 2002). In addition, efficient scattering 
is essential for the acceleration of cosmic rays. For instance, 
scattering of cosmic rays back into the shock i s a vital compo - 
nent of the first order Fermi acceleration (see Longair 1997). 
At the same time, stochastic acceleration by turbulence is en- 
tirely based on scattering. 

It is generally accepted that the energy of turbulence is most 
probably due to supernovae explosion^ and cascaded down to 
small scales, where resonance with CRs of moderate energies 
happerQ. 

For years magnetic turbulence has been treated in ad hoc 
manner, with variations which included plane Alfven waves 
moving parallel to magnetic field, a combination of these 
waves and so-called 2D perturbations, as isotropic magnetic 
turbulence with Kolmogorov spectrum and anisotropic turbu- 
lence with constant degree of anisotropy. Numerical simula- 
tion contradicted to all the models above, however. 

Unlike hydrodynamic turbulence, Alfvenic turbulence de- 
velops a scale-dependent anisotropic, with eddies el ongated 
along the magnetic field dGoldreich & Sridharlll995L hence- 
forth GS95). The GS95 model describes incompressible 
Alfvenic turbulence, which formally means that plasma p = 
Pgas/Pmag, the ratio of gas pressure to magnetic pressures is 
infinity. But it may be conjectured that GS95 scaling should 
be approximately true for m oderately com pressible plasma 
dLithwick & Goldreichl 1200 lh . Studies in (ICho & La zarian 
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120021 henceforth CL02) showed that the coupling of Alfvenic 
and compressible modes is weak and that the Alfvenic 
modes follow the GS95 spectrum in compressible medium. 
This is consistent with the analysis of observational dat a 
dLazarian & Pog osvan 2000; Stani mirovic & Lazarianl 2001). 

The turbulence injected on large scales may correspond 
to GS95 model and its extensions to compressible medium 
is less efficient in scattering of CRs compared to the esti- 
mates made assuming that magnetic turbulence consists of 
plane waves moving parallel to magnetic field. A cascade 
of Alfvenic perturbations initiated at large injectio n scales is 
shown to be really inefficien t for the CR scattering (IChandranl 
2000: 1 Yan & Lazarianll2002l) . 

At the same time, one should not disregard the possibil- 
ities of generation of additional perturbations by CR them- 
selves. For instance, the slab Alfvenic pert urbation can be 
created, e.g., v ia streaming instability (see WentzeJ Il974t 
ICesarsky |[l980). These perturbations are present for a range 
of CRs energies (e.g., < lOOGeV in interstellar medium for 
the streaming instability) owing to non-linear damping arising 
from ambient turbul ence (lYan & La zarian 20021 120041 hence- 
forth YL02, YL04, iFarmer & Goldreichl 120041) . Instabilities 
induced by anisotropic distribution o f CRs were also sug- 
geste d as a possibility to scatter CRs dLerchdl 19671 : iMelrosd 
IT971 . 

Recent work by Laza rian & Beresnyakl ((2006, henceforth 
LB06) proposed that compressible turbulence can induce gy- 
roresonance instability, which is an important feedback pro- 
cesses that can create slab modes to efficiently scatter CRs. 
This process is claimed to be more efficient in scattering CRs 
compared to direct action of turbulent fluctuations. 

This gyro-kinetic instability is induced by the anisotropic 
distribution of CRs, which is caused by the compression aris- 
ing from large scale turbulence. The degree of anisotropy, is 
determined by the compression on the scale of CR mean free 
path in LB06. In their treatment the growth of CRs needs to 
be balanced with steepening in order to prevent the diffusion 
from approaching the Bohm limit, where the mean free path 
of particles becomes comparable to the Larmor radius. In this 
paper, we shall investigate the nonlinear suppression of the in- 
stability by considering the self-adjustment of anisotropy due 
to scattering with the magnetic perturbations. The instabil- 
ity reaches a stabilized growth rate due to the feedback of the 
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increased perturbations on the anisotropy of CRs. 

In what follows, in §2 we give a brief background on tur- 
bulence and CR scattering, in §3 we introduce the gyroreso- 
nance instability of CRs, in §4 we formulate our approach to 
the problem. We outline our findings and the implications of 
our results in §5. In §6 we discuss our results. The summary 
is provided in §7. 

2. MHD TURBULENCE AND COSMIC RAY SCATTERING 

The propagation of cosmic rays (CRs) is mainly affected 
by their interaction with magnetic field. The properties of 
turbulence are vital for the correct description of CR propaga- 
tion. For instance, the scattering frequency is much dependent 
on the anisotropy of t urbulence dLerche & Sc hlickeise feOOlt 
lYan & LazarTanll2002l) . 

Numerical simulations (s ee ICho & Vishniad 1 20001 : 
Maron & Goldreichl l200lt iMuller & Biskamd 200o[ 
Cho. L azarian. & Vishniac 2002; ICho & Lazarianl 2002, 
20031 see also the book of iBiskamd 120031 as well as, 
Cho et al.ll2003tlEImegreen & Scaloll2004l for reviews), show 
that Alfvenic modes exhibit scal e-dependent anisotropy. This 
anisotropy was first discussed in Goldr eich & Sridhan dl995l 
hencefoth GS95). However, the understanding came later that 
this anisotropy is observable only in the system of reference 
related to the local direction of magnetic field on the scale 
of the eddy under con s ideration (|Laza rian & Vishniad[T999l: 
ICho & Vishnlacl l2000b iMaron & Goldreichl l2001l) ~~On the 
intuitive level it can be explained as the result of the following 
fact: it is easier to mix the magnetic field lines perpendicular 
to the direction of the magnetic field rather than to bend 
them. However, one cannot do mixing in the perpendicular 
direction to very small scales without affecting the parallel 
scales. This is probably the major difference between the 
adopted m odel of Alfvenic p erturbations and the Reduced 
MHD (see Bi eber et al.l 11994). The corresponding scaling 
can be easily obtained from the critical balance condition, 
i.e., kj^dvk ~ fen va, where feii and k± are the parallel and per- 
pendicular components of the wave vector k, va is the Alfven 
speed, Sv^. oc £j v is the turbulence velocity. Throughout the 
paper, || and _L are defin ed in reference t o the local magnetic 
field. As proven later bv lCho et al.l (120021) . the mixing motion 
perpendicular to magnetic field is essentially hydrodynamic 
and energy in the turbulence cascade is conserved locally 
in phase space. From these arguments, the scale dependent 

anisotropy fe|| °< k 2 ^ and oV^. oc k ± 1 ^ « fe,, can be derived. 

Owing to the existence of the parallel cascade, Alfvenic tur- 
bulence can cascade to the resonant scale of CRs, which is 
characterized by the parallel scale k^ res ~ £l/v, where Q. is 
the Larmor frequency and v is the particle speed. The scat- 
tering efficiency is, however, substantially reduced because 
of the anisotropy of the Alfvenic turbulence so th at we can 
completely neglect it (Yan & Lazarian 2002, 2004) except for 
ultra-high energy CRs. 

The distribution of energy between compressible and in- 
compressible modes depends, in general, on the way turbu- 
lence is driven. Naturally a more systematic study of different 
types of driving is required to determine the energy partition. 
In the absence of this, in what follows we assume that equal 
amounts of energy are transferred into fast and Alfven modes 
when driving is at large scales. 

While particular aspects of the GS95 model, e.g., the 
particular value o f the spectral index, ar e the subject of 
controversies (see Miiller & Biskamp 2000; Boldvrev 2006; 



iGogoberidzel 120071; iBeresnvak & Lazarianl 120091) . we think 
that, at present, GS95 model provides a good starting 
for developing models of CR scattering. MHD turbu- 
lence can be decompo s ed int o Alfven, slow and fast modes 
(ICho & Lazarianll2002l l200l . Slow modes are passive and 
follow the same GS95 scal i ng as do the Alfven modes 
dLithwick & Goldreichl 120011; ICho & Lazarianl 120021) . Fast 
modes, instead, are marginally affected by Alfven mo des and 
follow acoustic type cascade (ICho &L azarian 2003). They 
have been identified by YL02.YL04. lYan & Lazarianl (2008, 
hencefoth YL08) as the major source of CR scattering in in- 
terstellar and intracluster medium. 

Until recently, test particle approximation was assumed in 
most of earlier studies and no feedback of CRs is included 
apart from the streaming instability. Turbulence cascade is 
established from large scales and no feedback of CRs is in- 
cluded. This may not reflect the reality as we know the en- 
ergy of C Rs is comparab le to that in turbulence and magnetic 
field (see lKulsrudll2005l) . It was suggested by LB06 that the 
gyroresonance instability of CRs can drain energy from the 
large scale turbulence and cause instability on small scales 
by the turbulence compression induced anisotropy on CRs. 
And the wave generated on the scales, in turn, provides ad- 
ditional scattering to CRs. In this paper, we would like to 
provide quantitative studies based on the nonlinear theory of 
the growth of the instability to take into account more accu- 
rately the feedback of the growing waves on the distributions 
of CRs. 

3. GYRORESONANCE INSTABILITY OF COSMIC RAYS 

Plasma is known to be subject to numero us instabilities, on e 
of them being gyroresonance instability (iGary et al.l 119971) . 
LB06 suggested that an analogous instability operates in CRs 
and is important for CRs subject to compressions arising from 
MHD turbulence. 

3.1. Origin of growth 

In a nonuniform plasma, where magnetic field strength 
varies, instability can develop due to anisotropy of particles 
arising from the adiabatic invariant p, which is preserved 
when the variations happen on length-scale much larger than 
the particles' Larmor radii and timescale much longer than the 
Larmor period. 

To understand the essence of the instability it is important to 
recall that the definition for the magnetic moment of a charged 
particle in magnetic field is 



where q, p, ri are the charge, momentum and Larmor radius 
of the particle, c is the light speed, B is the strength of the 
magnetic field. For non-relativistic particles, p = t Ei C) ±/B is 
proportional to the kinetic energy of particles £^ j_; for rel- 
ativistic particles, however, p±v± deviates from the real ki- 
netic energy, p is an adiabatic invariant, indicating that p±v± 
changes with the strength of magnetic field. In the regions 
where magnetic field increases, the perpendicular pressure of 
the particles is higher than the parallel one, inducing mirror 
and gyroresonance instability; in the places where magnetic 
field decreases, the opposite is true, resulting in firehose in- 
stability. Waves are generated through the instabilities, en- 
hancing the scattering rates of the particles, their distribution 
will be relaxed to the state of marginal state of instability 
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even in the collisionless environment. We define the degree 
of anisotropy of particle distribution as 



where 



Wn 



W ± =J d 3 pf '(p)vxpj 
Wn = f d 3 pf(p)y m 



(2) 



(3) 



is equivalent to the perpendicular and parallel components of 
the energy density of CRs. 

In the presence of background compressible turbulence, the 
CR distribution is bound to be anisotropic because of the con- 
servation of the adiabatic invariant /j. Such anisotropic distri- 
bution is subjected to various instabilities. While the hydrody- 
namic instability requires certain threshold, the kinetic insta- 
bility can grow very fast with small deviations from isotropy. 
In this paper, we shall focus on the gyroresonance instability. 
The rate at which the anisotropy A varies is determined by 
turbulence compressional rate, 



1 dB 



5v ± 
va 



(4) 



where CO = v p h/l c is the frequency of the compressional modes 
at the compressional scale l c , 8v± °< f/. is the perpendicular 
component of the turbulence velocity. Since beyond mean 
free path, the anisotropy is effectively reduced because of 
scattering, only compression < X, the mean free path of CRs, 
can alter the degree of anisotropy. As 0)Sv ± /va « and the 
power law index of turbulence v is less than 1, the compres- 
sion rate in fact increases with the decrease of the scale. So 
different from LB06, we adopt the turbulence damping scale 
l c as the compression scale provided that l c ^> ri. For Kol- 
mogorov scaling, V = 1/3. In the case of X < l c , as in some 
collisionless medium, the compression is still supplied from 
l c , but with a reduced rate coSv^/va x ^ 1+v A,/'c- Using the 
compressional time CO -1 as the timescale for growth of the 
anisotropy, LB06 estimated the anisotropy A ~ Sv/va- We no- 
tice, however, that the growth time for the anisotropy should 
be constrained by the scattering time of CRs z scatt . If l scatt — 
X/c were adopted, the anisotropy of C Rs A ~ % xcatt dA /dt 
would not reach the threshold va/v (see Kulsrud 2005) un- 
less the mean free path of CRs is much larger than the turbu- 
lence compressional scale, which is limited by the damping^. 
In fact, this constraint does not apply because we are dealing 
with an evolving system rather than an equilibrium one. The 
wave amplitude e is growing at a rate F gr and so the mean 
free path is decreasing. Therefore the scattering rate is de- 
creased compared to the equilibrium value CI/ En = CIB 2 / 87ie. 
The reduced rate may be estimated as below. The variation of 
anisotropy A can be expressed as 



dA 
~dt 



vA = 



1 (dW\ W, dW 



Wn 



dt 



Wn dt 



' TgrEN/^CR, (5) 



where p^ff = 87tW||/B 2 is the ratio of CR pressure to mag- 

6 In collisionless envii 
of compressible modes 
IBrunetti & Lazarianll2T)0"7l) 



6 In collisionless environments, for instance, the dissipation scale 
of compressible modes can be very large I Yan & Lazarian] 120041 ; 



netic pressure. In the above equation, we utilized the fact that 
A << 1 and so W±/Wn ~ 1 and that during the isotropization 
induced by the growing waves, the loss in W± is greater than 
the gain in Wii so that wave can grow at the expense of the 
CRs' energy. From the above equation, the isotropization rate 
can be estimated as 



TgrEN 
PcflA' 



(6) 



By balancing the rate of decrease in anisotropy due to scatter- 
ing and the growth due to compression, we get 



pcgcoSv 

T gr VA 

X = r p /E N . 



(7) 



The magnetic field strength is altered by compressible 
modes. In low P medium, the compression is mainly from fast 
modes. For fast modes, their cascade is slower than their wave 
oscillation. This means that they have periodical oscillations 
unlike slow modes turbulence. A compression/expansion is 
followed by another, therefore a growth rate must be faster 
than the w ave frequency in order for the instability to grow 
(lYanll2009h . In high p medium, the situation is different and 
the contribution from slow modes is dominant. 

In the above equations, we did not account for collisions. 
This is physical because the scale we consider is less than the 
mean free path. In LB06, the only feedback was passively 
incorporated in X, which decreases as the wave grows. This 
process, however, will not stop until X reaches Larmor radius, 
namely, the Bohm diffusion regime. In this paper, we shall 
provide a quantitative treatment of the nonlinear suppression. 

3.2. Growth rate 

The distribution function o f CRs /(x,p,f ) is governed by 
the Vlasov equation (see, e.g. Kulsrud 2005): 



3/ 5/ 



E(x,f)- 



vxB 



dp 



0, 



(8) 



where v,p are velocity and the momentum of the particles 
respectively. E is the electric field. We can take a perturbation 
expansion for the distribution function 



/ = /0+/l+/2 + . 



fi+l < fx 



(9) 



if the field perturbation is small dB <C Bo, where Bq is taken 
in z directioifl 

Assuming a monochromic Alfven wave propagating along 
z direction, the linear perturba tion f\ can be obtained using 
the quasi-linear theory (see e.g. Kulsrud 2005). The lineariza- 
tion of the Vlasov equation is 



- ico/i + ikvji + - (v x B ) ■ V p /i = -q 



By integrating it, one can get 



E 



v x (k x E) 



CO 



•V, 



(10) 



7 The gyroradii are much smaller than the turbulence injection scale unless 
one is considering ultra-high energy CRs. So the condition SB <g Bo is well 
justified for moderate energy CRs. 



4 



Yan & Lazarian 



iqv±E^ ( i dfp 



CO 



<°- k \\ v \\ d/o 
V_L 9p± 



(11) 

where /j + corresponds to the right-moving cosmic rays and 
/j~ is for left-moving cosmic rays. CO, k are the wave frequency 
and the wave number. £2 = £2o/Yl is the relativistic Larmor 
frequency, £2o = qB/m/c, m is the proton mass. From the 
Maxwell equations, one can obtain 



Ei 1 



k 2 c 2 
CO 2 



471/ J i 
CO 



= 0. 



where 



Hi J d 3 v\fi 



(12) 



(13) 



is the current density arising from the perturbed distribution 
f\. Then, the dispersion relation can be obtained through 
eas. dl 1112113b 



k 2 c 2 



CO 



(14) 



and the perturbed part of dielectric tensor 
47t/Ji 2%q 2 



K 



± _ ' 



coEi 



co z 



d 3 p 



co — fen vii ±Cl 



kv±- h (co — feiiVl 

dpu " 

(15) 

where the subscript 'j' refers to different species. For CO <S £2, 
K° = c 2 /v\, representing MHD waves parallel to the mag- 
netic field with coo = kvA- Insert CO = COo + COi into the above 
equation, one can find the growth rate corresponding to the 
imaginary part of the perturbed Kj, 



= 3(©i) 



coovj 
" 2c 2 



3(*f) 



(16) 



The growth rate of the gyroresonance instability can be eas- 
ily derived from this equation: 



2 2 

■■% e vA 



CO 



ticrip > mQ.a/k\\) 



{A-A cr )Q, 



f(A-A cr )(l 



l-a 



(17) 



where n is the gas number density, co p = y 'Anne 2 Jm is the 
proton plasma frequency, n cr °< p l ~ a is the number density of 
cosmic rays, and (see e.g., LB06) 

Acr =Ya/v, 

e = ^(a + 2)(a-l)Bf(3/2,a/2), 
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4 x 10- 10 cm- 3 



n C r{r P > ro) 



pc, (18) 



in which ro is the g yro-radius of lGeV particles, B t represents 
a beta function (see Abramowitz & Stegun 1965). 

The sign of Y gr determines whether the wave grows ('+') 
or damps ('-'). It depends on the bracketed expression in 




FIG. 1. — Anisotropy of particle distributions and waves excited by them. 
The oblate distribution on the left can be caused by rarefaction/expansion of 
the field and the prolate on the right can be realized due to compression of 
the field. The dashed lines represent the resonant momenta p res . The dashdot 
circles are t he paths along which the derivative of f needs to be considered 
(see Eq |19l . As explained in the text, for each distribution, right moving 
and left moving CRs excite right-handed and left-handed polarized (annular 
arrows) waves moving in opposite directions (the lines with arrows). 



dp± 



Eq.dT5t. which is a de rivative of / a long a circle centered at 
= (seelulsrud 20Ql, 

dp± co — kv z 
dpu kv± 



(19) 



The wave grows if there are more and more particles along the 
way it propagates. Unlike with shifted distribution (streaming 
instability), only circularly polarized waves can be excited by 
the gyroresonance instability. First we decide the directions of 
wave vector resonant with right and left moving cosmic rays, 
respectively. Then we determine the polarization of the waves 
bearing in mind that CR protons always rotate clockwise with 
respect to the magnetic field. For an oblate distribution of par- 
ticles with p± > pn, the left moving, right circularly polarized 
wave can resonate with the right moving CRs (with super- 
script '+') and gain energy from them as it encounters more 
particles in the wave propagation direction (see Fig^ight). 
So do the right moving, left circularly polarized waves with 
the left moving CRs (with superscript '-'). For a prolate dis- 
tribution, the opposite is true. They are illustrated by FigQ] 

4. FEEDBACK PROCESSES AND STEADY STATE GROWTH 
SOLUTIONS 

4.1. Nonlinear suppression 

In this section we shall provide a rigorous study of the sup- 
pression of the growth of the instability. The wave-particle in- 
teraction was t reated using the second order theory (SOT) in 
iGarv & Tokarl (11985b for bi-Maxwellian distribution of ther- 
mal particles. These particles have Maxwellian distribution 
characterized by Zj|,7j_, two different temperatures parallel 
and perpendicular to the local magnetic field. We shall mod- 
ify the SOT for the CRs' power-law distribution, i.e. f(p) « 
p- a - 2 . 

Assuming small-amplitude fluctuation, we can easily ex- 
tend the Vlasov equation to the second order. The slowly 
varying second order distribution of particles due to the 
stochastic interactions with the electro-magnetic perturba- 
tions can be described by the ensemble average of the second 
order extension of the Vlasov equation, 
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d <f 2 > gvxB d <f 2 > 



dt 

= -e< (Ei + 



c 

v x Bi 



8p 
8p 



IT 3/0 



>, 



(20) 



The evolution of the anisotropy of the CRs A can be derived 
from this equation. The electric field Ei is perpendicular to 
the magnetic field and therefore increases the perpendicular 
energy of CRs. The magnetic perturbation Bi scatter parti- 
cles and thus only causes exchange of particle energies in the 
parallel and perpendicular direction. The second moment of 
Eq.d2"0b gives 



and 



dW_i(k) 
dt 

ggjlW 
dt 



2(a + 2) 



■■s±(k) 



de(k) 
dt • 



(26) 



(a+l)(a- 
(a+2) 



-co 



pr 



-2k 2 c 2 



S\\{k)=2 



~dT 

dWu 



2q<Ei- / \fid 3 p > 



2q 



<B, 



at c 



d pv± x v/i >, 



d i py± x v/i >, 
(21) 



a+ l)(a- 
a(a- 1) 

(a + 2; 
(oH 



-&Cr{p > Pres) + 1 



-co„ 



2) 



-fc 2 c 2 



The two terms on the right hand side of the first equation cor- 
respond to the acceleration and scattering. They are given by 
(see Appendix for the detailed derivation), 

q<Ei-J vfd 3 p>= J ^|J|2 3 [® 2 ®*(K+e- +K~e+)] , 

| < B t • y d 3 p\ ± x v/i >= -/ d*™ {<»pr (r + e 

+ 3 [(co,- + n)co 2 /T + e_ + (co,- - i2)co 2 /r e+] } , 

where G> pr = ^/4-Kn cr (p > p res )q 2 /m is equivalent of plasma 
frequency, but for CRs, co, is the real part of CO. 



_(a+l)(a+3) 
From this, we can get 

dt \Wi 



/{\c»\ 2 + c 2 k 2 }, 

Pc«(p > Pres) + 1 



/dt- 



1 

wii 



[S±-(A- 



er e 



(27) 



(28) 



By equating Eqs.(l4ll. ([28l . we get 



re 4 



■l)Sf gr £ N = 0, (29) 



e± = 



E* ± E ± 

871 



For a symmetric distrib ution, w e expect that e 4 
r + = r~. Then from Eqs. (12 11221 ). we can get 



(23) 
£_ and 



Since SinSj. are of order 1 (see Fig.l), we see that the results 
(22)fr° m above equation is comparable to our earlier estimate in 
Eq.(|7]i- There are two unknowns A,En in the equation d29l ). 
There has to be another relation to determine them. We as- 
sume here that the system relaxes to the marginal state of in- 
stability, and A = A cr = va/v. For the compression by fast 
modes, co = va/I c -, §v± ~ V(l c /L) v , V the injection velocity 
of turbulence at the injection scale L. We get for the compres- 
sion by fast modes, 



dW±_ 
dW\\ 

~~dT 



dkS ± (k) 
dkSu (k) 



dejk) 
dt ' 
de(k) 
dt ' 



e N = $ C rM a 

I [5,|(l+A)-5j , ifX>l c 



U_ 

l c A 



a-l 



(30) 



where 



e(*) = 



8tc 



k 2 c 2 

= ( 1 +l^ir)( e + + e -) 



(24) 



(25) 



£n = < &crM A j- 



l c A 



a-l 



/ [5|,(l+A)-5j} 2 , ifX<l c 



(31) 



is the total energy density of the wave and S±,Su are mea- 
sures of the scattering efficiency, determining how fast is the 
response of the distribution of CRs to the wave perturbations. 
They increase with the decrease of wave frequency (see Ap- 
pendix for the detailed derivation). In general, their function 
is to reduce the anisotropy. 

For gyroresonance instability, the growth of the wave hap- 
pens at the gyroresonance frequency of the CRs, k res . There- 
fore in the first order approximation the above evolution equa- 
tion can be differentiated, and accordingly we obtain (see Ap- 
pendix for detailed derivation) 



where Ma = V / va is the Alfvenic Mach number, v — 1/3 for 
the Kolmogorov turbulence and v = 1/4 for the Iroshinikov- 
Kraichnan turbulence. We adopt Kolmogorov scaling^] for the 
fast modes in this paper. 

The parallel damping scale of slow modes in high p 
medium is the proton diffusion scale $ l / 2 l m f D , where l m f „ 
is the mean free path of thermal particles (Ba rnes! Il966l) . 
For slow modes in high p medium, CO = va/I c , 5vj_/va = 
Sv/v^cosO ~ (I c /Ia)~ x I 2 , where Ia = L/M\ in the case of 
Ma > 1. Inserting it into Eq.d2"9l>, we get 



8 At a given scale I c , the compression and wave growth would be stronger 
if IK scaling is used. 
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FIG. 2. — Scattering coefficients Sj_,5ii vs. CR energy. 



1/AfbM 



1/r„ 



1/r„ 



I [S||(1+A)-Sjl] , ifk>l e (32) 



eN= {* CR k{t) u^i 

I [S,|(l+A)-5 ± ]}2 , //A</ c (33) 

4.2. bottle-neck for the growth due to energy constraint 

The creation of the slab waves through the CR resonant in- 
stability is another channel to drain the energy of large scale 
turbulence. This process, on one hand, can damp the turbu- 
lence. On the other hand, it means that the growth rate is 
limited by the turbulence cascade. The energy growth rate 
r^e cannot be larger than the turbulence energy cascading 
rate, which is 1 /2pV 4 /va/L for fast modes in low p medium 
and pv^/U for slow modes in high p medium. This places a 
constraint on the growth: 



M*v A /L, p<l, 



va/Ia, 



P>1, 



Thus the upper limit of wave energy is given by 



M 2 A Li/(LA)(r p /ro) a -\ p<l 
Z./MXAv/ro) 01 - 1 , p>l. 



(34) 



(35) 



The growth is induced by the compression at scales < A. 
Therefore, in the case that r gr EN reaches the energy cascad- 
ing rate, fast modes are damped at the corresponding maxi- 
mum turbulence pumping scale Xft, = r p /z^. If Xft, is larger 
than the original damping scale l c , then l c should be replaced 
by X/b, and one needs to iterate until the results are self- 
consistent. 

4.3. Summary of Procedures 

The gyroresonance instability can be driven by CRs with 
anisotropic distribution which is induced by compressible tur- 
bulence. Fast and slow modes dominate the compression in 



FIG. 3. — The spectral energy density of slab waves that is transferred from 
the large scale compressible turbulence via the gyroresonance instability of 
CRs. In the case that the instability grows up to the maximum energy rate al- 
lowed by the turbulence cascade, large scale turbulence is truncated at Xf/, and 
the wave amplitude E(k)dk ~ e'^ is given by Eq. )35t . Note that the picture 
is different from LB06, namely, the feedback on the large scale turbulence 
occurs only in some cases when the scattering is not sufficient to prevent the 
waves from growing to the maximum values. 

low and high p medium respectively. The compressional rate 
is maximized at the smallest scale-damping scale of the tur- 
bulence. Since CRs do not 'remember' the perturbation of 
magnetic field and its anisotropy beyond its mean free path, 
compressions larger than the mean free path of CRs will not 
induce anisotropy. Therefore, we adopt min(A, l c ) as the com- 
pressional scale that drives the anisotropy of CRs. We first 
calculate the damping scale of the compressible modes. If the 
damping scale is larg er than the original mean free path Ao, 
one should use the Eq. (1301321 ) to obtain the wave amplitude 
and the new mean free path A = r p jz^. If the dam ping sc ale 
is smaller than the mean free path Ao, we use Eqs.( l3 1133b to 
get 8jv and A = r p jz^. 
The resulting energy growth rate r ?r e should be compared 

with the large scale turbulence energy cascading rate V 3 /L. 
If it is larger than V 3 /L, the growth is limited by the energy 
budget, rather than the nonlinear feedback. Then eq.d35ll can 
be used to estimate the growth (see Figf3]l. 

The growth rate should also overcome the damping by 
background turbulence (YL0 2, 04; Farmer & Goldreich 2004, 
iBeresnyak & La zarian 2008). If the growth rate is lower than 
the damping rate 



r f / = Vim I \jLMfp, 



(36) 



the instability does not operate. Lm,Vlm are the injection 
scale of strong MHD turbulence. Lm = Ia, Vlm = Va if 
M A > 1, and L M = LM 2 A ,V LM = v A M 2 A if M A < 1. 

For fast modes, one also needs to take into account the tur- 
bulence compression rate at the dissipation scale co = Va/Ic- 
If the growth rate is smaller than the compression, the subse- 
quent expansion will smear out the anisotropy, the trigger of 
the kinetic instability, and the growth will not happen. 

We show results of our study in figures H] [5] |6] In the 
plots, we also demonstrate several relevant dimensionless pa- 
rameters: degree of anisotropy A and the ratio of CR pres- 
sure to magnetic pressure, Pck, as well as the growth and co. 
The mean free path of CRs are plotted together with the ear- 
lier estimates from LB06. Compared to LB06, the degree of 
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anisotropy A is much smaller, and set by the marginal state 
of instability. This results in a much lower growth rate, and 
therefore smaller range of energies of CRs that can induce 
the instability because of the constraint by damping by back- 
ground turbulence. In addition, the nonlinear feedback due to 
the scattering by the growing waves is more efficient in sat- 
urating the growth at small amplitude (e# < 1) compared to 
the feedback considered in LB06. We do not need additional 
wave steepening to stop the wave from growing to large am- 
plitude (or the diffusion to reach the Bohm limit). Our treat- 
ment is self-consistent. 

5. ASTROPHYSICAL IMPLICATIONS 
5.1. Importance of the instability 

Our present work confirms the finding in LB06 that the gy- 
roresonance instability plays an important role for the prop- 
agation of cosmic rays in turbulent astrophysical fluids. We 
note, that the coupling of the compressible and incompress- 
ible motions reported in Cho & Lazarian (2003) between 
compressible and Alfvenic modes is strong enough at the 
scale of driving for diverting about 20% of energy into com- 
pressible cascade, even if the driving is purely solenoidal. 
While we still do not know many details of the MHD cas- 
cade (see Beresnyak & Lazarian 2009), the fact that a notice- 
able portion of energy is going into large scale compressions 
is generally accepted. It is those compressions that the gy- 
roresonance instability requires to induce substantial changes 
of the dynamics of the astrophysical fluids. 

In view of the above understandings, it looks wrong that 
present day numerical codes ignore the effects of cosmic rays. 
We may expect an appreciable change of the compressible 
dynamics if cosmic rays drain energy from the compressible 
motions. The consequences of the process for the ISM and 
molecular cloud formation will be discussed elsewhere. 

At the same time, the plane parallel Alfven waves with 

3> are expected to interact efficiently with cosmic rays^J 
Apart from scattering this should induce cosmic ray accelera- 
tion via second order Fermi process. The consequence of this 
important process is still to be evaluated quantitatively for dif- 
ferent astrophysical environments. 

Our quantitative study shows that LB06 in their simpli- 
fied treatment overestimated the growth rate and therefore the 
range of energies of cosmic rays over which the instability is 
active. The reason is that the degree of anisotropy A is over- 
estimated in LB06. The nonlinear feedback that we study in 
this paper limit the growth through the scattering and there- 
fore A is maintained at the level of marginal instability. Due 
to the same scattering, our steady state wave energy/mean free 
path is smaller/larger than that in LB06 (see the right panels 
of FigJU [5] |6}. The quenching of growth they considered by 
only reducing the mean free path is not adequate, and in fact 
does not apply here as we consider the compression at the 
minimum scale of turbulence instead of the mean free path. 
Turbulence energy cascading rate is another constraint that 
sometimes is more stringent than the scattering FT. 

5.2. Particular cases of astrophysical environments 

9 In the Solar and Magnetospheric community it is accepted to talk about 
the energetic particles rather than cosmic rays. This distinction is not impor- 
tant for the basic processes we discuss. 

10 The effect on turbulence cascade was considered in LB06, but not the 
feedback on the wave growth. 



5.2.1. Galaxy 

In interstellar medium, compressible turbulence is sub- 
jected to various damping processes. In the case of ionized 
phase, the damping of fast modes varies with the wave pitch 
angle 9 (YL02, 04, 08). The turbulence truncation scale in- 
creases with the pitch angle. Approaching to 90°, the trun- 
cation scale is the largest partially because of the field line 
wandering (YL04, 08, see figHfe/f), where the mean free path 
of cosmic rays X can be smaller than the truncation scale l c . In 
the mean time, at smaller pitch angles, there is less damping 
and the mean free path X is larger than the truncation scale. 
In the low p medium, the fast modes are the dominant source 
for the compression of magnetic fields and we focus on fast 
modes for the low (3 Galactic halo. Using the compression 
by the fast modes at the mean free path, we obtain a fairly 
large growth of the waves £o with amplitude larger than that 
allowed by the energy budget of the turbulence (see Fig0k). 
Turbulence is thus damped at the maximum mean free path of 
CRs Xmax = r p/z l N- Compression below this scale is in fact 
reduced and eqj30|should be adopted to recalculate the ampli- 
tude of the wave. We show the results in Figj4] 

Another case we consider is hot ionized medium (HIM), 
which differs from halo in the sense P > 1. In the high 
P medium, the dominant compression originates from slow 
modes. In the collisionless medium, the damping scale of 
slow modes, the proton diffusion scale is quite large. We use 
therefore Eq.d33t to calculate the wave amplitude. We check 
the consistency by comparing the final mean free path with 
the damping scale. The results are shown in Fig [5] 

5.2.2. Clusters of galaxies 

Observations of galaxy clusters show nonthermal compo- 
nent together with the thermal component in the intracluster 
medium. Turbulence may be induced by cluster of mergers 
and accretion of matter at the virial radius. Both observation 
and simulations suggest an appreciable amount of energy re- 
sides in turbulence. The typical velocity of turbulence at the 
injection scale is expected to be around 500-1000 km/s, which 
is highly super-Alfvenic because of high plasma p. In this 
case, the magnetic field correlation length Ia could be smaller 
than the coulomb mean free path. The effective mean free 
path should be the minimum of the two if not accounting for 
further scattering of thermal particles due to plasma instabili- 
ties. 

Because of the similarity with the HIM, the results for intr- 
acluster medium (ICM) can be obtained with the similar pro- 
cedure. They are shown in Figs|6] Note that because of the 
plasma instabilities, the real thermal proton mean free path is 
less than the one we use. This means that the mean free paths 
for CRs we obtain is an upper limit and practically, they can 
be even smaller. 

6. DISCUSSION 

We demonstrate that the steady state amplitude of the slab 
wave increases with the CR energy. This is because among all 
the parameters, the growth rate changes y l L ~ a ) most rapidly 
with the energy. And since the scattering rate linearly depends 
on the growth rate, the wave amplitude, inversely proportional 
to the scattering rate, increases with the energy of CRs (eqs|6] 
IT]). In the case that the growth is constrained by the energy 
budget of turbulence, the wave amplitude has a simple relation 
e iV 06 Tl 1 x k l ~ a on the energy (see Eq[35]). 
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T(K) 


n(cm 3 ) 


B(juG) 


M A 


L(pc) 


n cr (10- lu cm- 3 ) 


l c (cm) 


X fb (cm) 


halo 


2x 10 6 


O.OOl 


5 


1 


30 


4 


2xl0 19 


1.84 x 10 17 


HIM 


10 6 


0.005 


2.5 


1 


30 


4 


4.9 x 1() 16 




ICM 


10 8 


0.001 


1 


10 


3 x 10 5 


4 


4.2 x 10 19 





TABLE 1 

The physical parameters we adopted for different medium. HIM=hot ionized medium, ICM=general intracluster medium. 
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FIG. 4. — (a): the amplitude of the growing wave using different approaches. ,o is obtained with pumping at the mean free path, which exceeds the maximum 
value eJJ, allowed by the energy of turbulence and is therefore unrealistic; E r j is calculated using the reduced compression from the damping scale set by the 
maximum mean free path obtained through ,o; the real amplitude is the minimum of Eai = min(e](,,£„/) (see text for the details), (b): The ratio of cosmic ray to 
magnetic pressure (5„ , the degree of anisotropy A and the wave energy of the slab modes at the steady state vs kinetic energy of CRs in Galactic halo, (c): the 
steady state wave growth rate, turbulence compressional rate and damping rate. Growth only happens if the growth rate is larger than the damping rate. For the 
compression by fast modes, compressional rate is another threshold, (d): the mean free paths, obtained in this paper and in LB06. 
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We show that the gyroresonance instability operates only 
for CRs of a limited energy range. On the low energy end, 
the turbulence energy can limit the growth of the instability 
(as in the case of halo). On the high energy end, the growth is 
constrained by the feedback of the growing slab waves on the 
anisotropy of the CRs' distribution. The mechanism we dis- 
cussed earlier (YL02,04,08), namely, CR scattering by large 
scale fast modes are, therefore, still dominant for higher en- 
ergy CRs. 

Indeed the damping by background turbulence is the ma- 
jor constraint for the growth of the instability. In the envi- 
ronments where Alfvenic turbulence is weak or damped, the 
instability can extend to higher energies. For instance, in the 
regions where the ionization degree is low, the slab wave aris- 
ing from the instability will protrude to a scale of tj^ /L l l 3 
(LB06). 

In the calculation on collisionless plasma, we used simply 
the coulomb mean free path. The effective viscosity in colli- 
sionless plasma should be less because of additional scattering 
of thermal particles caused by plasma instabilities. The con- 
sideration of this complex processes is beyond the scope of 
this paper. However, one thing we can be certain is that with 
the reduced viscosity, higher compressional efficiency can be 
achieved and the mean free path of CRs is further reduced and 
CRs are better confined, which is crucial for the acceleration 
of CRs. 

Our results show that turbulence compression provides an 
additional channel for CR scattering and the feedback of CRs 
is important for the cascade of turbulence, confirming LB06's 
findings although our quantitative conclusions are very much 
different. In future MHD simulations, it is necessary to take 
into account the contribution from CRs. 

7. SUMMARY 



We have considered the gyroresonance instability of CRs 
in MHD turbulence. It arises from the compression from 
the compressible turbulence. The growth of the instability 
is self-adjusted due to the feedback of the slab wave on the 
anisotropy of CRs. We find that 

1. Compression by large scale turbulence results in 
anisotropic distribution of CRs, which induces gyrores- 
onance instability. 

2. The growth of the instability is balanced by the feed- 
back of the slab waves on the anisotropy of CRs. The 
high energy cutoff is determined by the balance be- 
tween the growth rate and the damping of the wave by 
the large scale turbulence. 

3. The slab waves generated through gyroresonance in- 
stability limits the degree of the anisotropy of CRs, 
and is an efficient isotropization mechanism for CRs< 
lOOGeV. 

4. The feedback of the instability on the large scale tur- 
bulence should be taken into account. In some cases, 
turbulence cascade is truncated because of the drain of 
energy to the small scale waves. 
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APPENDIX 

DERIVATION OF SCATTERING EFFICIENCIES S± ,S 

Define 

co - kvg d/b dfo 



D= s^L+k-^. (Al) 

v_l dpx_ dpu 
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2e3i < E ■ T > 
-2e3i<E- j d 3 p\fi > 

Size 2 



-e T 3 / d*p 



co fevii-co=Fi2 



-D 



r 3(co 2 co*/r ± )e: 



(A2) 



where we applied Eg. ( 111115b and the fact < R < A(x,t)B(x,t) >= 1 /29t[A*(k,f)fi(k,f)]. We use K to express the result as it is 
directly associated with the growth rate (EqJTSl 



2e 
c 



3i< J dVr[vxBi]/i > 



-2e<R < / d 3 p\ ± 



-/V x E 



co 



/i > 



e5R < / d 3 pkvi 



v -E 



co 



-/i > 



= -e T 3 / t/ 3 p 



87te 2 Dvi^vn 
|co| 2 fevii -CO^H 



-s 3 / rf 3 ^M ( ! 



j3 STte 2 3/, 



co±ft 



fevii - COT^ 



= -e T 3 / d i p-^v ± ^3((a) + ri x3[(j) 2 ((j)±Q.)K ± ]e. 



|co| 2 3p 



CO 



(A3) 



For a power law distribution fo °< p a 2 , vj_dfo/dpj_ = — (<x + 2)cp 2 L p a 5 . Insert it into the 2nd to the last term into Eq. (IA3b . 
we get 



3 87te 2 3/o 
16ne 2 



=F (a+l)(a + 3)|co| 2 

2(a + 2) 
(a+l)(a + 3) 



P|| \ 4n cr{P > Pres) 



3 (co) 



P|| \ 4co pr (p > 

kol 2 



3(co)e^ 



Therefore, 



-— < Bi • / c/ J pv ± x v/i > 



(A4) 



(A5) 



From Eqs. d21IA2IA6t , we get 

3W ± 



2(a + 2) 
(a+l)(a + 3) 

2(a + 2) 



Pres 



,2 



±1 4e T 



co 



CO 



dt 



-4£ 
4E 



(a+l)(a + 3) 
2(a + 2) 



(a+l)(a + 3) 



CO' 

R 

,2 



(A6) 



where £ is performed over the right moving CRs ('+') and left moving cosmic rays ('-'). In the above equations, 



or 



: (co, + iy) 2 = (co 2 + y 2 ) + 2/co r r, 



K ± =c 2 /v 2 (l-2ir ± /& ). 
Taking into account that y <C CO,- and co r » coo, we get 



(A7) 
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W±(k) 
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where we used the fact that de/dt = 2T gr (l +A: 2 c 2 /|co| 2 )(e+ +e 
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